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Abstract 

For a quantum system governed by a non-Hermitian Hamiltonian, we studied the 
problem of obtaining an optimum Hamiltonian that generates nonunitary transfor- 
mations of a given initial state into a certain final state in the smallest time r. 
The analysis is based on the relationship between the states of the two-dimensional 
subspace of the Hilbert space spanned by the initial and final states and the points 
of the two-dimensional complex Bloch sphere. 
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In view of recent results on optimal quantum evolution and its possible rela- 
tion with quantum computation and quantum information processing, there is 
increasing interest in the quantum brachistochrone problem(see e.g., jH, 0,0,0] 
and references therein). The problem consists of finding the shortest time r 
to evolve a given initial state {ipi) into a certain final state {ipf) under a given 
set of constraints. Using the variational principle, Carlini et al [l| has shown 
that a Hermitian Hamiltonian r has a nonzero lower bound. Later, Bender et 
al j3| demonstrated that, for non-Hermitian PT-symmetric quantum systems, 
the answer is quite different, and the evolution time r can be made arbitrary 
small, despite the fact that the eigenvalue constraint is held fixed as it is for 
the corresponding Hermitian system. 

Recently, Brody and Hook have considered the same problem as in [ij 
using the symmetry properties of the quantum state space and without em- 
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ploying a variational calculus. The purpose of this paper is to extend the 
approach of Brody and Hook to quantum systems governed by non-Hermitian 
Hamiltonians (for details and discussion of non-Hermitian physics see, e.g., 
(sl, @, 0, @, @, Eo] ) • We explore the non-Hermitian time-optimal evolution prob- 
lem from the geometric viewpoint, using the relationship between the states 
of two-dimensional subspace of the Hilbert space 7i, spanned by and 
\ipf)i and the points of two-dimensional complex sphere (the so-called com- 
plex Bloch sphere). 

In the quantum brachistochrone problem, finding the shortest possible time 
requires only the solution of a two-dimensional problem. Namely, it requires 
finding the optimal evolution time for the quantum system governed by the 
effective Hamiltonian acting in the subspace spanned by and \il)f) [21,0, Si- 
Let the set {I'i/'i), |'?/'o), (V'il) (V'ol} form the bi-orthonormal basis of the two- 
dimensional subspace spanned by the states \ipi) and I?/'/): 

= (^olV-o) = 1, and (V^.lV'o) = (V^ol^^) = 0. (1) 
Then the final state can be written as 

\iljf) = a\^i) + b\iJo), {^f \ = d{ij,\ + b{i^o\, (2) 

and a generic non-Hermitian Hamiltonian H acting in the two-dimensional 
invariant subspace of the Hilbert space Ti. reads 

H =X,{\tp,){ij,\ + |^o)(V^o|) + ^cos9\[\ij,){iJi\ - \iJo){H 

+ ^sme (e-^'^l^,) (^ol + e^lV'o) , (3) 

Ao being a complex parameter and = A+ — A_, where A+, A_ are complex 
eigenvalues of the Hamiltonian H. Thus, the time-optimal evolution problem 
can be formulated as follows: determine the complex angles 6, (p for which, 
under a given set of eigenvalue constraints, the initial state iV'j) evolves into 
the final state \il)f) in the smallest time r. 

Let \il){t)) = ui{t)\tlji) +U2{t)\t/jo) and {^{t)\ = ui{t){iji\ + U2(t)'(^o| be a 
solution of the Schrodinger equation and its adjoint equation 

z^^im) = H\m), = i^(t)\H (4) 

where we have chosen the units with h = 1. As can be easily shown, the vectors 
\uit)) = uiit)\u^) + U2{t)\ui) {u{t)\=u,it){u^\+U2it){ui\, (5) 

where = ( ^ j and = ( ° j denote the up/down states, respectively. 
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satisfy the following Schrodinger equation: 

= Hef\u), -^g-^{u{t)\ = {n{t)\H,f. (6) 
The effective Hamiltonian entering in ([6]) is given by 



H, 



'a„ \ 




V A„j 





COS 6 e *'^sin6' 
e**^ sin 6 — cos 6 



(7) 



and can be expressed in terms of the Pauli matrices as 



H,f = Xol + ^n-(T, (8) 

where 1 denotes the identity operator, and f2 = (fi sin 6 cos (p, Q sin 6 sin ip, Q cos ( 
Using complex Cartesian coordinates, 

X = Q sin 6 cos ip, Y = Q sin 6 sin ip, Z = Q cos 6, X,Y, Z E C, (9) 

where Q, = ^/X^ + Y'^ + Z^, we obtain 



(10) 



It is then easy to show that the complex Bloch vector defined as n = {u\(t\u), 
where cr are the Pauli matrices, satisfies the complex Bloch equation 

^ = O X n, (11) 




which is equivalent to the Schrodinger equation ([6]) (see, e.g., [ll|, ll2|). The 
vector n(t) = {ni(t) , n2(t) , n^it)) , being a complex unit vector, traces out a 
trajectory on the complex 2-dimensional sphere S"^. In the explicit form, its 
components in terms of the states ([5]) read 

ni = U1U2 + U2U1, 712 = i{uiU2 - U2U1), 713 = UlUl - U2U2. (12) 

Hence, we obtain the Bloch vector corresponding to the initial state l^/^j) as 
= (0,0,1). 

Let Wf == {uf\cr\uf) be the complex Bloch vector associated with the final 
state l^pf), and let r be the amount of time required to evolve the initial 
state \ipi) into the final state {ipf). Then, using the time dependent solution of 
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^ sin 6 cos ^^(1 — cos ^t) cos (p + sin 6 sin Q,t sin ip ^ 



n(t) 



sin 6 cos ^^(1 — cos Qt) sin ip — sin 6* sin Qt cos ^9 
cos^ 9 + sin^ cos Qt 



(13) 



with rij = n(0) = (0,0, 1) and rif = n(r) 
find 



(sin X cos 7, sin X sin 7, cos x) , we 



-'^ cot el tan - -iWtan^^-tan^^ 



e*^ cot ^ tan - + iWtan^ ^ - tan^ ^ 



cos 



cos^ I — cos^ 6 



sin 6' 



sm 



sin I 
sin ^ 



This yields 



cos((y9 — 7) = cot tan 



X 



tan ■ 



fir 



sm 



cos^ I — cos^ 9 



(14) 
(15) 

(16) 

(17) 
(18) 



and we obtain the evolution time as r = where 
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arctan 



sm 



/ cos^ I — cos^ 6 



(19) 



In addition, since r is a real positive function, the following constraint must 
be imposed: arg\E' = 0. This implies 



arg Q = arg arctan 



sin I 



' cos^ I — cos^ 9 



(20) 



We note that Eq. fl20|) can be considered in different ways. For instance, for 
a given Q, it restricts all possible final states to those that have the angle 
X satisfying this equation. In contrast, if the eigenvalue constraint is given 
by = const, then fl2Ul) can be considered as an equation to determine the 
argument of Q. Finally, for a given Q and x, the equation (120!) yields an implicit 
function dependent on the variables ^9 and '^9, and it should be considered 
in finding the evolution time. 

For a given Q, using (fT8l) and (120|) . we find that the critical points of the 
function t{9) are defined by the following equations: 
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(21) 
(22) 



Then, from Eq. (!2T]) . we obtain two solutions: ^6 = 7r/2 (536' = 0) and '^6 — > 
±oo. Inserting = 7r/2 into Eq. (!22l) . we find argfi = argx, and, for the 
second solution, we have ^6 = 7i/2 ± [aigQ — argsin(x/2)). 



Further study of the critical points shows that there is no solution with a finite 
value of \6\ yielding the minimum of the evolution time. Indeed, when argf2 = 
argx, the first solution, 6 = 7r/2, is related to the saddle point for the function 
r = t{9), and the second one yields r ^ while '^9 ±oo (see Fig. [1]). Thus, 
for a quantum mechanical system governed by a non-Hermitian Hamiltonian, 
the evolution time r has a zero lower bound and may be taken to be arbitrarily 
small. Moreover, since, for any finite value of |^|, the minimum of r does not 
exist, the non-Hermitian Hamiltonian cannot be optimized. However, for a 
quantum mechanical system governed by a Hermitian Hamiltonian, the latter 
can be optimized. Indeed, in this case, we have = 0, and, hence, the 
saddle point becomes the point of a local minimum (Figs. [HE]). This agrees 
completely with the results obtained in [H, [ij. 




Fig. 1. The evolution time r is depicted as a line of intersection of the sur- 
faces z = \^\ and arg^* = 0. Left panel: O = 1, x = tt. Right panel: 
Q, = 1 + 0.25i, X = -|- 0.25i. The eigenvalue constraint is given by O = const. 



Applying (jUl) to eliminate (p from the effective Hamiltonian (171), we obtain 
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Fig. 2. Plot of evolution time r vs. ?R.9 and Q6. Left panel: = 1, Qil. = 0.1, 
X = TT + i. Right panel: = IjSri = 5,x = tt. The constraint is taken as 
= const. 

i^e/ = Hef{^, 9, X, l)- In the explicit form, we have 

QcosO ( 1 e-^T(tanf -iJtan^^-tan^f) 

i/e/ = AqI H , , , 

^ \^ e'^ ( tan f + zy'tan^ 9 - tan^ | J -1 

(23) 



To obtain the Hamiltonian in terms of initial and final states, we need to solve 
Eq. ([2]) and determine lipo). After a lengthy calculation using the explicit 
solution of the Schrodinger equation (jl]), we obtain: 

cos — - i cos 9 sin — j e"*^"", 6 = -ie-'^°^e"^ sin ^ sin — , (24) 

cos — + zcos^sin — e*^"^, 6 = ze^^o^e"*'^ sin ^ sin — . (25 

2 2 / ' 2 ^ ^ 

Then, applying ([12]), ([H]) - ([IH]) and ([MD, (ESD, we obtain 

ii" =Aol + „ . ^. 2 y ( fcosgJcos^^ -cos^6' + isin^)|?/'/)(^i| 
2 sm b' sm f \ ^ ' ^ 

+ (^cos^-^cos2- -cos2^-zsin-j|^/'i)(V'/| j - 2T (^1^^) + 

(26) 

where cos^ | = (^^l?/'/). 

The generic optimization problem on the complex Bloch sphere is, for a given 
final vector n/ = (sin x cos 7, sin x sin 7, cos x) , to determine the complex angle 
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6 for which, under a given set of eigenvalue constraints, the initial vector 
rij = (0, 0, 1) evolves into the final vector Uj in the smallest time r. As follows 
from the preceding analysis, the non-Hermitian Hamiltonian can be optimized 
only if '^O = 0. In particular, it emerges that the Hermitian Hamiltonian is a 
special case of the latter when, in addition, $57 = '^x = 0- 

As can be observed in Figs. [H [2l the point ^ = 7r/2 is a saddle point for 
the evolution time r = T{x,y), where x = ^6 and y = Q6. For 6 = 7r/2, 
the straightforward computation yields r = and, in addition, we have 

argfi = argx. This relationship implies that, for a given Q, only final states 
rij satisfying argx = aigQ can be reached from the initial state rij. Inserting 
6 = tt/2 into (123!) . we obtain the effective Hamiltonian as 



To compare with Hermitian quantum evolution, let us note that, for Hermi- 
tian quantum systems, all variables (x, d, f) and Q are real. In that case, the 
point 6' = 7r/2 is the point of the local minimum, and the Hamiltonian fl28l) . in- 
deed, becomes the optimal Hamiltonian. This is in accordance with the results 
obtained in the previous works of Carlini et al and Brody and Hook [l], 0] . 

In what follows, we consider some illustrative examples, starting with the 
generalization of the results obtained by Carlini et al in their paper [l| on the 
quantum brachistochrone problem. In jlj, the Hamiltonian constraint has been 
imposed on the standard deviation AH of the Hamiltonian. In the standard 
normalized state, we have (Aif)^ = {ip\H'^\'il)) — {{ip\H\ip)Y = const. The 
natural generalization of the complex energy variance to the non-Hermitian 
Hamiltonian is given by (AE)^ = {'il)\H'^\%l)) — ((■?/' | if |-?/;))^, where {i>\ip) = 1. 
Then, applying this to the Hamiltonian ([S]) and with the help of Eq. (1131) . we 
obtain AE = (l]sin^)/2. 

Imposing the Hamiltonian constraint as AE = const, from Eq. f|T8|) we obtain 




(27) 



Finally, from Eq. (l26l) . the "optimal" Hamiltonian is found to be 




(28) 



tAE 



sin I 



(29) 



tan 



sin 6 
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Fig. 3. The evolution time r is depicted as a line of intersection of the sur- 
faces z = \^\ and arg^ = 0. Left panel: 0, = 1, x = 7r/2- Right panel: 
Q = 1 + 0.25i, X = 7i"/2 + 0.25i. The eigenvalue constraint is given by AE = 1. 



This yields r = |\^|, where 

sin 6^ / sin I \ , ^ 

^ = — — arctan ^ = (30) 

\^cos2f -cos^e; 

In addition, since r is a real function, one should impose the following con- 
straint arg \^ = 0. In Fig. [3], the evolution time r is depicted as the line of 
intersection of the surfaces z = and arg\^ = 0. 

The following example is related to the recent controversy on the possibility 
of achieving faster evolution in a quantum mechanical system governed by 
a non-Hermitian PT-symmetric Hamiltonian as compared to the equivalent 
Hermitian system [i, 13, 14, 15 1. Note that the critique of the results obtained 
in Ref. [3] is essentially based on the following theorem 15 1: The lower bound 
on the travel time (upper bound on the speed) of unitary evolutions is a univer- 
sal quantity, independent of whether the evolution is generated by a Hermitian 
or non-Hermitian Hamiltonian. Analyzing the proof of this theorem, one can 
see that it is based on following assumption: the minimal travel time is real- 
ized by quantum evolution along the geodesic path in the Hilbert space joining 
initial and final states. This is true in the case of the Hermitian Hamiltonian; 
however, as we further show, for a non-Hermitian Hamiltonian, the situation 
is quite different. 

Without loss of generality, we can further confine our attention to the case of 
= 0- Taking into consideration that eigenvalues of PT-symmetric Hamilto- 
nian are real, and employing Eq. ([TB|1 . we find ^9 = tt/2 and 3?((/9— 7) = ±7r/2. 
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In what follows, to be definite, we will choose the upper sign. Then, denoting 
the imaginary part of 6 as rj and substituting 6 = tt/2 + irj into Eq. ([8]), we 
obtain 

1 



H. 



Anl + - f2 ■ cr 

" 2 



(31) 



where O = f2(cosh?7cos<y9, cosh?7sin<^, — -isinhr^). Further, we assume f2 > 0. 
Using a convenient parametrization. 



X 



VL cosh 7] cos (p, y = Q cosh rj sin ip, z = Q sinh rj 



(32) 



we have x"^ + y"^ — z'^ = . Thus, the eigenvalue constraint Q = const defines 
a one-sheeted hyperboloid embedded in a fiat three-dimensional space M^. 
Finally, inserting 9 = 7i/2 + irj into Eq. (fT3|) . we obtain 



^ -i sinh r] cosh 77(1 — cos ^It) cos ip + cosh rj sin fit sin y9 ^ 



n(t) 



— i sinh ?7 cosh — cos fit) sin ip> — cosh r/ sin Qt cos 
cosh^ ?7 cos ilt — sinh^ r/ 



(33) 



and, evidently, n(0) = (0,0, 1). 



We show further that, for a given Q, the evolution of a quantum mechanical 
system governed by the non-Hermitian Hamiltonian with real spectra can be 
considered as motion over the surface of the unit one-sheeted hyperboloid 
EP. Let n{t) be the Bloch vector given by (l33ll . Then, we define the map 
n(t) m(t) as 



m 



m 
m 



rt" = cosh^ 7] cos fit — sinh^ r] (34) 
ra"*^ sin — ra^ cos = cosh 7] sin fit (35) 
i{n^ cos {p + sin ip) = sinh r] cosh //(I — cos fit) (36) 



This yields gijm^m^ = 1, where Qij = diag(l, 1, — 1) is the indefinite metric, 
and, hence, Eqs. f lM|) - fl36|) define the map 5^ H^. Thus, the vector m(t) 
traces out a trajectory on the unit one-sheeted hyperboloid H^, while the 
quantum mechanical system evolves in the Hilbert space. From flM|l - fl36|) . it 
then follows that the initial Bloch vector iij = (0, 0, 1) maps to rrij = (1, 0, 0), 
and, for the final state defined by nj = (sin x cos 7, sin x sin 7, cos x), we have 

Uf ^ nif = (^cosx, \l sin^ x + tanh^ 77(1 — cosx)^, tanhr7(l — cosx)) (37) 



The amount of time r required to evolve the initial state rrij into the final state 
m/ can be found from Eq. (ITSl) by substituting 6 = 7i/2 + iri. The computation 
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yields 



n 



arctan 



sm 



cos^ ^ + sinh^ r] 



(3J 



Therefore, it follows that a non-Hermitian Hamiltonian with real spectra r has 
the upper bound t < where x/^ is the evolution time for the optimal 

Hermitian Hamiltonian (!28l) and can be obtained from (138!) taking 77 = 
(see Fig. H]). Our results are in agreement with those obtained previously by 
Bender et a/ in [3], and they contradict the conclusions on the impossibility 
of achieving a faster evolution for non-Hermitian quantum evolution than for 
Hermitian evolution 15 . 




Fig. 4. Plot of evolution time r vs. x V = foi' ^ non-Hermitian Hamiltonian 
with real spectra (0, = 1). As can be observed, for a given x, the evolution time has 
a maximum at the point r] = 0. 



To compare r with the geodesic evolution time r^, we must consider the in- 
duced intrinsic metric ds"^ = Qij dm^dm^ on the unit one-sheeted hyperboloid 
H^, where gij = diag(l, 1, — 1). Let m be a unit vector with respect to the 
indefinite metric on Qij. Using a convenient parametrization, 

= cosh p cos u, = cosh p sin u, m? = sinh p (39) 

where —00 < p < 00 and < u < 2it, we find 

ds"^ = cosh^ p di/ - dp^ (40) 
Any geodesic on can be obtained as a solution of the geodesic equation 
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cPx^ ■ dx^ dx^ , , 

+ rU-^-^ = o (41) 



ds^ ■'^ ds ds 
For the metric fl4UI) . it takes the form 



d'^p . . f dh'\^ , , 

+ smhpcoshp — =0 (42) 
ds^ \ds J 

d^p , dpdu , , 

d? + 2tanhp^- = (43) 

We note that, setting m = (cosh p cos z/, cosh p sin z/, sinhp), one can rewrite 
the geodesic equation as 



^ + m = (44) 

Its solution is given by 



duii 

m(s) = cos s rrij + sin s — — (45) 

ds 

Taking the initial state as rrij = (1,0,0) and the final state as uif from Eq. 
(1371), we have 



m = (^ cos s, y 1 + tanh r^tanh — sin s, tanh?7 tanh — sin sj (46) 

Now, comparing fllBl) with Eqs. flM]) - fl5B]) . we conclude that m(t) determined 
by the Schrodinger equation yields the geodesic evolution only if 77 = 0. In 
addition, for the geodesic motion, we have s = fit. Using these results, we can 
rewrite (l38l) as 

2 , sinf 
T = Tg— arctan ^ (47) 

X J cos^ I + sinh^ rj 

where Tg = x/fl is the geodesic evolution time. This imphes that r < if 77 7^ 
0, and T = Tg if 7] = 0. Furthermore, whereas, for 6 = ^/2 {'^6 = rj = 0), the 
optimal Hermitian Hamiltonian provides the shortest travel time, for a non- 
Hermitian Hamiltonian, the geodesic motion over the one-sheeted hyperboloid 
yields only the upper bound on the evolution time. 

Computing the distance on between the initial nij and final lUf states, 
calculated along the curve C defined by Eqs. ( l34l) - ( !36l) . we obtain 

L = f \l gijih^rh^ dt = tQ cosh 77 = 2 cosh rj arctan I ^ j (48) 

-'^ \ J cos^ I + sinh^ 77 / 
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As can be easily shown, L is bounded as 2sin(x/2) < L < x-, where the 
upper bound is determined by the geodesic distance Lg = flTg = x between 
the initial and final states on (Fig. ED- The lower bound, L 2sin(x/2), 
is reached while — > oo. Note that the geodesic distance between nij and 
lUf calculated on on with the help of the indefinite metric PUI) coincides 
with the geodesic distance between the related initial n, and final states on 
the conventional Bloch sphere S'^. However, in the latter case, the system is 
governed by the optimal Hermitian Hamiltonian, and = x is the shortest 
distance on S"^ between rij and rij. 




-5.0 



Fig. 5. Plot of L vs. X ^iid rj for the quantum evolution generated by the non-Her- 
mitian Hamiltonian (|31|) . 



To compare the evolution time for the non-Hermitian and Hermitian Hamil- 
tonians and conclude which one is faster, one should not only fix the geodesic 



distance between initial and final states, as was pointed out in jl5|. but also 
impose the same set of constraints for both cases. Here, following [ill, we as- 
sume that the Hamiltonian constraint 2AE = Q cosh t] = const is held. Then, 
denoting ficosh?] as u, we obtain from Eq. (1471) the following result: 

2 rising 
T = — arctan , ^ (49) 

^ ^uo^ - sin^ I 

This yields 

- sin - < r < - (50) 

As can be seen, the evolution time r reaches it maximum Tmax = x/^ when 
Q = u. This implies t] = 0, and we have the quantum mechanical evolution 
governed by the optimal Hermitian Hamiltonian for which 6 = 7t/2. The lower 
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bound, Tmin = (2/c<j) sin(x/2), is obtained at the exceptional point defined by 
the condition f2 = 0. Similar consideration of the distance (l48l) between nij 
and nif yields 



2uj sin f 

uiT = — — arctan , (51) 
^ Juj^ - sin^ f 



and 



2sin|<L<x (52) 



The upper bound Lmax = being identical to the geodesic distance between 
the initial and final states defined either on the Bloch sphere S"^ or one-sheeted 
hyperboloid H^, is achieved for the Hermitian Hamiltonian. The lower bound, 
Lmin = 2sin(x/2), is obtained at the exceptional point. 

Thus, as has been shown above, the geodesic distance Lg between rij and 
Uf calculated over the Bloch sphere S"^ is identical to the geodesic distance 
between rrij and mj computed over the one-sheeted hyperboloid H^. Moreover, 
the amount of time Tg required to evolve iij into rij on the Bloch sphere and 
that required to evolve rrij into on by the geodesic evolution, is the 



same. This is in accordance with the conclusions made in 15 . However, in the 



case of the Hermitian Hamiltonian, Tg is the lower bound on the evolution time, 
and, for the non-Hermitian Hamiltonian, it yields only the upper bound on the 
evolution time. Furthermore, as follows from Eq. (1481) for the non-Hermitian 
Hamiltonian, the evolution speed v = ds/dt is given hj v = u = Q cosh rj (we 
recall that t] = Q9). Hence, v > Vg, where Vg = fl is the speed of the geodesic 
evolution. Similar consideration of the quantum mechanical system governed 
by the Hermitian Hamiltonian yields v = QsinO, and, obviously, v <Vg. This 
proves that, indeed, non-Hermitian quantum mechanics can be faster than 
Hermitian evolution 



In summary, we have formulated a geometric problem on the complex Bloch 
sphere to find the optimal non-Hermitian Hamiltonian and the optimal time 
evolution for a given pair of initial and final states and eigenvalue constraints. 
In contrast to the Hermitian quantum system, generic non-Hermitian Hamil- 
tonians generate non-unitary transformations {ipj) such that the evo- 

lution may be realized in an arbitrarily short time. This is in agreement with 
the previous results obtained for PT-symmetric quantum systems and some 



specific non-Hermitian Hamiltonians [4, [14 . 
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